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Abstract 

The Af = 2* Snper-Yang-Mills theory (SYM*) nndergoes an inhnite se- 
qnence of large-A qnantnm phase transitions. We compnte expectation 
valnes of Wilson loops in fc-symmetric and antisymmetric representations 
of the SU{N) gange gronp in this theory and show that the same phe¬ 
nomenon that canses the phase transitions at hnite conpling leads to a 
non-analytic dependence of Wilson loops on k/N when the conpling is 
strictly inhnite, thns making the higher-representation Wilson loops ideal 
holographic probes of the non-trivial phase strnctnre of SYM*. 


1 Introduction 

The path integral of any Af = 2 gange theory on S^ can be calcnlated P by 
means of the snpersymmetric localization [2]. Localization reduces the path 
integral to a matrix model giving us direct access to strong coupling, and in 

*Also at ITEP, Moscow, Russia 
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particular to the regime of interest for holography, when N is inhnite and 
the’t Hooft coupling A = is large. 

The SYM* theory is particularly well-suited for this purpose, since its 
holographic dual is explicitly known |3], while the Af = 2* localization matrix 
model can be solved at strong coupling in El E] by more or less standard 
methods of random matrix theory [7j. Various predictions of holography can 
then be confronted with ab initio evaluation of the held-theory path integral 
m Perfect agreement found so far has left, nevertheless, one feature un¬ 
explained from the holographic perspective. Localization predicts an inhnite 
sequence of quantum phase transitions that occur at certain critical values of 
the’t Hooft coupling in the large-V limit PIID]. Large-iV phase transitions 
are of course very common in matrix models m , but here the real critical be¬ 
havior occurs such that correlation length goes to inhnity. The ensuing phase 
transitions also have £nite-iV counterparts [12]. Similar transitions occur in 
other localization path integrals uniiis], and it would be very interesting to 
understand their holographic origin. 

Unfortunately, tuning the’t Hooft coupling to some given critical value 
is very difficult in holography. Well-established holographic methods work 
at strictly inhnite coupling, leaving us little hope to observe the transitions 
as discontinuities of correlation functions in A. However, in the SYM* the 
n-th critical coupling grows with n as Ac"’^ 2 ; uni, and most of the phase 
transitions do happen when the coupling is very large. Certain observables 
may then have criticality-induced non-analyticities even when the coupling 
is strictly inhnite. Identifying such observables is the goal of this paper. 

The phase transitions are driven by nearly massless states that give rise to 
singularities in the eigenvalue density of the localization matrix model. The 
eigenvalue density can be directly probed by D-branes in the dual supergrav¬ 
ity background na. The classical D-branes see the coarse-grained density, 
in which the resonances are smoothened out [9], but D-brane quantum huc- 
tuations may be sensitive to the non-trivial phase structure of SYM*. 

Another type of observables calculable by localization are Wilson loops. 
They might appear useless for detecting phase transitions at hrst sight. There 
are two scales in the problem, M and IjR, where M is the mass scale of the 
SYM* theory and R is the radius of compactihcation on 5^. For the most 
part of this paper we regard the radius of as an IR regulator and thus 
assume that MR » 1. At strong coupling a new scale p ~ \/XM » M 
emerges. Wilson loops in the fundamental representation probe the smallest 
scales of order l/R, the D-branes probe the largest scale of order p, while 
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resonances that cause critical behavior occur at the intermediate scale of 
order M. In this paper we consider Wilson loops in higher representations 
of the gauge group, and show that by varying the size of representation it is 
possible to scan the whole spectrum of scales, including the resonance region. 

Wilson loops of rank k ^ N are dual to D-brane conhgurations carrying 
the fundamental string charge [151 El EH ESI [la EHl EU, and have been 
extensively studied in the W = 4 super-Yang-Mills theory [Ml EH EH EH EH > 
for which the matrix model is just Gaussian [25l EH]- far as M = 2 theories 
are concerned, higher-rank Wilson loops have been calculated for the Af = 2 
superconformal QCD [27]) whose solution at strong coupling is known [2H| 
but is very different from that of SYM*. 


2 Wilson loops 

The SYM* theory has the same held content as AA = 4 SYM - gauge helds, 
six adjoint scalars and four Majorana fermions, but its Lagrangian includes 
explicit mass terms for half of the helds (four scalars and their Af = 2 super¬ 
partners). We denote this common mass by M. 

The Wilson loop in representation TZ is dehned as 


Wn{C) 


P exp 



ds {ix^A^ + |i:|<I)) |, 


( 2 . 1 ) 


where $ is a scalar from the vector multiplet, coupling to which makes the 
Wilson loop locally supersymmetric and well-dehned in the UV. 

Localization computes the Wilson loop for the equatorial contour of 
[T]- The partition function on the four-sphere of radius R localizes to an 
eigenvalue integral: 


Z= f ^aY[ Zi~ioopiai-aj) e (2.2) 

i<j 


The integration variables are the Coulomb moduli that parameterize space- 
time independent zero modes of the scalar that enters the loop operator: 

<I)o = diag(ai,...,aAr). (2.3) 
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The non-zero modes of <h, as well as all other degrees of freedom have been 
integrated out, leaving behind a one-loop contribution [Ij: 


-2^1-ioop( 2 ^) ~ 




H{x + M)H{x-M) ’ 


00 / 

n=l ^ 


1 H -;;— j e " 




(2.4) 

The Wilson loop expectation value is obtained by simply replacing quan¬ 
tum fields in (2.1) by their values on the localization locus and subsequently 
averaging over the Coulomb moduli: 

Wn{C) = . (2.5) 

We denote by L the length of the contour C. The average is now dehned by 
the eigenvalue partition function (2.2). 

When C is the big circle on 5"^, this result is exact (up to instanton 
corrections neglected here because of the large-iV suppression), in which case 


L = 2ttR. 


( 2 . 6 ) 


Localization, strictly speaking, is not applicable to other contours, but the 
leading exponential behavior for asymptotically large Wilson loops should be 
universal and largely insensitive to the contour’s shape. We thus expect that 
our results apply to any sufficiently big contour, including contours on IT^, as 
long as ML » 1. This assertion is based on universality and has less rigorous 
grounds compared to localization, but has been checked holographically for 
Wilson loops in the fundamental representation |3]. 

We concentrate on ra uk -fc symmetric or antisymmetric representations: 


-fc (2.7) 



and will eventually consider the scaling limit iV ^ 00 , fc ^ cx), in which the 
ratio 

/4 P.8) 

is kept fixed. Such Wilson loops are dual to macroscopic D3 or D5 branes, 
depending on whether representation is symmetric or antisymmetric 
nziiiniEo], that carry k units of electric flux on their world-volume. On the 
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field-theory side, the fc-symmetric/antisymmetric Wilson loops are described 
by statistical mechanics of free bosons/fermions, for which the matrix eigen¬ 
values play the role of energy levels |23]. Let us review how this statistical 
interpretation arises in the matrix model. 

The fc-symmetric or antisymmetric characters are conveniently packaged 
into generating functions: 


= X! e 

k=0 


■‘"trK.e'", 


(2.9) 


For the antisymmetric representations, the sum terminates at k = N. An 
SU{N) character is a class function, in other words, a symmetric function 
of the eigenvalues Sn of the matrix E. When expressed in terms of eigen¬ 
values, the generating functions can be explicitly written as Bose and Fermi 
distribution: 

N 


( 2 . 10 ) 


n=l 


where (minus) the eigenvalue plays the role of the energy level and -u plays 
the role of chemical potential. 

Using this representation, we hnd for the Wilson loop (2.5) in the k- 
symmetric/antisymmetric representation: 


= W^ = L 


b 


C+iiT 

I 


dv 

27ri 


, kLu 


n[i 


^ e 


L{aj-u) 


r 


( 2 . 11 ) 


For symmetric representations, the contour of integration should be chosen 
such that C > aj for any j. For antisymmetric representations, C is arbitrary. 
We have rescaled u hj L for future convenience. 

The vacuum distribution of eigenvalues aj is characterized by the eigen¬ 
value density: 



( 2 . 12 ) 


At large-A^ the density does not fluctuate, because the action in (2.2) is 
0{N'^) while there are only N variables of integration. The weight in the 
integral representation for the rank-/c Wilson loop (2.11) is exponentially 
large, but the exponent is 0{N) rather than and the insertion of the 

Wilson loop does not affect the saddle point distribution of eigenvalues to 
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the leading order in 1/N. The statistical average in (2.11) is then replaced 
by the average with respect to the hxed ensemble determined by the master 
held p(a:)|3 


C+ilT 


f ^ e 

J 27rt 


NLT^{v) 


(2.13) 


where 


C-iTT 


r- 


-(u) = fu ^ - J dxp{x)\n{l^ , (2.14) 




and /, dehned in (2.8), is assumed to be of order one in the large-iV limit. 


The contour integral in (2.13) is of the saddle-point type, and to the 
leading order in l/N we hnd: 


\nW^ = + O {N^) , 

where u is determined by minimizing the free energy: 


f 


-M 


dx p(x) 

, L{u-x) ^ 


= /• 


(2.15) 


(2.16) 


This is the standard relation between canonical and grand canonical ensem¬ 
bles in statistical mechanics, with -v playing the role of chemical potential, / 
the density of particles and p{x) the level density. As expected, the exponent 
in the expectation value of the Wilson loop is proportional to the length of 
the contour. To compute the coefficient of proportionality, we hrst need to 


solve for the eigenvalue density of the matrix model (2.2), then hnd v from 


(2.16) and substitute the result into (2.14). 


The eigenvalue density for the localization matrix model is not known in 
general, but at large R and large A, a systematic perturbative solution can be 
constructed. The leading-order strong coupling solution was known for long 
time from the D-brane probe analysis of the dual supergravity background 
and can be easily reproduced from the saddle-point equations of the 
matrix model |1]. It has the form of the Wigner distribution: 


Poo(x) = — -\fp: 


2 


TT/i^ 


(2.17) 


^This is not quite true for the symmetric representations in certain range of parameters, 
as we shall discuss in later. 
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Figure 1: The eigenvalue density for /r = 7.3M. 


with the width proportional to the square root of the’t Hooft coupling 


/i 


s/\M 

2ti 


( 2 . 18 ) 


The leading-order is way too simple to describe the critical behavior of 
the model, which arises because of the complicated short-distance structure 
of the density, that at the leading order gets averaged over. Interestingly, the 
first strong-coupling correction fully reveals the short-distance singularities 
that are responsible for quantum phase transitions at finite A [SIE]. 

Although we expect from holography that the strong-coupling expansion 
goes in powers of 1/VA, the hrst correction to the density appears at the 
relative order ~ The correction has a very irregular, spiky 

^This results holds in the infinite volume. The finite-i? expression for the partition 
function on is obtained by replacing M with \J + l/i?^ g]. 
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structure (fig. 0 EE]: 


p{x) 2 ; 


TT/r^ 






M 


ttV 2/i^ 
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where {•} denotes the fractional part. The zeta-function (C(l/2,z) has a 
singularity at zero, and this produces an array of infinite spikes at resonance 
points X = ±/i T nM. 

The physical origin of the spiky structure in the density is due to the res¬ 
onance on massless hypermultiplets. The number of resonances that £t into 
the interval (-/r,/i) depends on the ratio 2p/M and each time 2p crosses an 
integer multiple of M, a new pair of resonances appears, leading to a fourth- 
order quantum phase transition pun]. The transitions are sharp only in infi¬ 
nite volume, at hnite R the spikes are rounded up and the transitions become 
smooth crossovers. The expression above assumes that the infinite volume 
limit is taken prior to the strong-coupling limit, in the sense that MR » \/\. 
If the strong-coupling limit is taken first, the spikes get damped, albeit at 
a very slow rate, and completely disappear distance M‘^R ~ pMRjy^ away 
from the endpoints [S]. 

In the bulk of the eigenvalue distribution, the comb-like structure of reso¬ 
nances is a small correction, which moreover disappears upon averaging over 
a sufficiently wide interval because 


1 

/<'«C(^.5)=0, (2.20) 

0 


However, near the endpoints the spikes are no longer suppressed compared 
to the averaged density. Indeed the two terms in (2.19) become comparable 
for p - X ~ M. The whole expression is actually not applicable near the 
endpoints. The edge behavior of the eigenvalue density is described instead 
by [5] 

^ [ ^—x ] 

1 



(p-x- 1), 


( 2 . 21 ) 


where [■] denotes the integer part. This function, shown in fig. smoothly 
matches with (2.19) at p - x » M. 















pM 



Figure 2: The density near the endpoint of the eigenvalue distribution. 


Extremely close to the endpoint, at fi-x ^ IjR, the structure of the peak 
starts to be resolved. The density has the following form in this regime [5]: 


p{x) = 


2M 





R 


dn 




K 


COS 2tik sin^ tik T^ (k) e K-i)-2nR(,.-x)K^ (^2.22) 


and is shown in Eg. This regime exists only for the theory defined on 
S^. The density near the edgepoints is known even without assuming that 
MR » 1 [5], but the general expression is substantially more complicated. 

With the density at hand, we can now compute the Wilson loop expec¬ 
tation value from (2.16), (2.15), (2.14). Depending on the parameter /, the 


main contribution to the integrals in (2.16), (2.14) will come from different 


range of x’s. We may identify three distinct regimes: 


Bulk regime: p-x ^ 0(p). In this case, the average density (2.17) is 


a good approximation. The spikes constitute a A correction as clear 


from (2.19). 


• Endpoint regime: fi-x ^ 0{M). The spiky structure of the density 
then appears in the leading order (|2.21 ). 


9 

























20 


40 


60 


80 


100 


2nR(ii-y.) 


Figure 3: The finite-i? resolution of the first peak, 
peak unresolved is shown in the purple line. 


The density (2.21) with the 


Non-universal regime: fi - x ~ 0(l/i?). The spike singularity is 
resolved according to (2.22). This regime is only dehned on as the 
density depends explicitly on the radius of the sphere, even though MR 
is still assumed to be large. 


The hrst two regimes are universal, in the sense that the results should apply 
to any sufficiently big contour and do not depend on the compactihcation 
(the radius drops out from the expressions for the density). The last regime 
requires that the theory is compactihed on the sphere, and the computations 
in that case only apply to the Wilson loop running along the big circle of S^. 


3 Antisymmetric representations 


The Wilson loop in the antisymmetric representation is expressed through 
the Fermi distribution, for which 1/L plays the role of temperature. When 
L is large, the effective temperature is low. Then (2.16) and (2.14) can be 
simplihed: 


J dxp{x) = / 


(3.1) 
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and 


dxp(x)x. (3.2) 

U 

This approximation is jnstified when the scale of variation of the density 
is mnch larger than 1/L, which is true in the bnlk and endpoint regimes, 
becanse pL » 1 and ML » 1, bnt not in the non-nniversal regime when the 
density varies on scales of order 1/i? = 271 jL. 


3.1 Bulk regime 

We begin with the bnlk regime, when the density is the same as in the 
Ganssian matrix model. In that case the answer is the same as in W = 4 
SYM np to rescaling of A by a factor of ML/27r, and it reads (THl 123] 


NML 


InW- = 


yX . 3 

- <51T' 


sim^ 


(3.3) 


where 9, defined as cos 6* = z///i, is the solntion of a transcendental eqnation 

1 


9 — sin 29 = tt/. 


(3.4) 


The correction term to the bnlk density in (2.19) is of order (9(A"^/^), 
bnt contribntes to the Wilson loop vev at relative order A"^/^, becanse the 
small-scale variations of the density average to zero. We are not going to 
stndy this contribntion in detail, since it is snbleading to the effect prodnced 
by the so far nnknown order A~^/^ correction to the density. 

From (3.4) we see that the bnlk regime corresponds to / = 0(1). When 


/ becomes smalQ the Fermi level v approaches the edge of the eigenvalne 
distribntion and eventnally the endpoint regime sets in. 


3.2 Endpoint regime 

Let ns define dimensionless variables u = {p - x)lM and v = {p- z/)/M. It is 


convenient to express the endpoint distribntion (2.21) for different intervals 


^or very close to one - all the formulas are invariant under / 
symmetry of antisymmetric representations. 


1-/ charge conjugation 
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of M, i.e. 


p{u) = ^ 


l 2M 1 
Try \/u 
1 /2M 


if M e (0,1] 

it “E (1.2] 


(3.5) 


-\[^Y!k-o /, ^ if M e + 1] 

I^ttV sjk+u-n 'i ’ J 


It is then clear that / should scale as {Mj ~ A To make this scaling 
manifest we introduce 

/ = -( — )'/ = 4v/^A-t/, (3.6) 

71 \ p / 

and assume that / - 0(1). 

The integral in (3.1) can be easily performed by summing over the con¬ 
tribution of different intervals of u that are smaller than v. For example, for 
n - 1 < n ^ n, with integer n, we have: 


1 ^ ^ 
^/n(^^) = y' p{u)du + - + J p{u) 

0 n-1 


du 


(3.7) 


where we have labeled / with a subindex n as a reminder of the domain of 
V. For the rescaled variables, the explicit forms are: 


/l(^) 

= 

(3.8) 

hiv) 

= s/v + Vn - 1 

(3.9) 

Mv) 

= Vn - 2 + \/n - 3, 

(3.10) 


and so on. The function /(n) is continuous but discontinuous in its hrst 
derivative at the critical points, which are its values at integer points, i.e. 

/'critical,n ~ /'n(^)- 


/critical,n - jl, 1 + ^^2, 1 + ^/2 + \/ 3 , . 


(3.11) 


Fig.i shows a plot of this function for 0 < n ^ 3. 


The free energy at the saddle point is given by (3.2). We are interested 


in expressing it in terms of the free parameter /, therefore, hrst, we must 
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Figure 4: The rescaled function f{v) plotted. The dash lines indicate the critical 
values /critical,n for n = 1, 2, 3. 


invert the function f{v). The inverse functions of /i(n) and f 2 iv) admits a 
simple closed form: 


0<n^l: Vi{f) = p 

(3.12) 

(p + l)^ 

l<n^2: V2if) = ~ 

4/2 

(3.13) 


V 3 {f) can be also expressed in radicals, but the expression is lengthy, and we 
do not display it here. 

Now, in terms of v{f), and using the saddle point equation (3.1) reex¬ 
pressed in the endpoint variables u and n, the free energy as a function of / 
is: 

^(/) 

(3.14) 


T (f) = - M‘^ J du p{u) u 


At strong coupling the hrst term is dominant, and to the hrst approximation 
the free energy is simply T~ = pf = pkjN. The Wilson loop in the rank-A; 
representation therefore behaves as which is just the /c-th power of the 

fundamental representation, as can be expected for small representations on 
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account of the large-A^ factorization. The second term is of relative order 
O(Mlij) = 0(l/'/X) and in the dual holographic description can be inter¬ 
preted as a quantum correction on the D-brane worldvolume. This term has 
singularities at the critical points (3.11), which means that the worldvolume 
effective held theory may have a non-trivial phase structure even at inhnitely 
strong coupling. 

Therefore, the phase structures are smooth, due to the subleading term, 
see hg. 

The phase transitions are of the second order, as can be seen by examining 
the derivatives of free energy. Using that dfjdv = Mp{v), we hnd: 


V 

The inverse density is singular for integer values of v, as can be seen in 
fig- HI and takes a hnite value to the left of the critical point going to zero to 
the right. The second derivative of the free energy consequently experiences 
a hnite jump across the phase transition. We can see this explicitly by 
computing the free energy on the hrst two intervals: 


dj^- 


= /i - Mu, 




P( 




ttM 

3A^ 


X s 


4/3 

P + 6 / 


0</^l 

l<f^l + V2 


(3.15) 


And it is straightforward to check that T is continuous at / = 1 together 
with its hrst derivative, while the second derivative experiences a hnite jump. 


3.3 Non-universal regime 

Extremely close to the endpoints, the hrst peak is resolved, and we need to 
solve (2.16) with the density from (2.22). From (2.16) we can see that the 
Wilson loop enters this regime at 

/ " 


i?l/2^3/2 


(Mi?)-fi2y-3/4^ 


which is an additional factor pMR smaller compared to (3.6). 
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f 


Figure 5; The subleading term of the free energy plotted, rescaled according to 
(3.15). The dash lines indicate the critical values /critical,n for n = 1,2,3. 


Considering again symmetric and antisymmetric representations in par¬ 
allel, we can see that the free energy scales as 




\/^M 


7r4(Mi?)3/2A3/4 




The scaling fnnction h^{i) can be written in the parametric form: 

oo 

h^(f) = -f IngT r ^ cos27rfi:sinVfi:r2(K)e-2'^(^'^'^-^) 
J K.2 


|ln(lTg) 


1 + K 


2 T’i( 1 , 1 + K]2 + K]±q) 


/ dhh , Q T— ,9 / \ 

—= COSZTTKSm TIKL [Kje' 

Jk 


■2Mi(ln K-1) 


X--2T'i(1, 1 + k; 2 + k; ±g). 

1 + K 


(3.16) 


(3.17) 


(3.18) 


The parameter q is the fngacity variable related to the chemical potential 
V in (2.16), (2.14) by g = For the antisymmetric representations 
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(fermions), q changes from zero (low density) at to infinity (high density), 
while for bosons the chemical potential must be negative and so q is always 
smaller than one. 

The low density approximation, q ^ 0, corresponds to Boltzmann statis¬ 
tics, when the difference between symmetric and antisymmetric representa¬ 
tion disappears. The hypergeometric function in the above equations can 
then be replaced by 1. We thus get: 

f = Ti^e~^q (<?-^ 0), (3.19) 


and 

3 

h^(f) = f (l-lng) = f In-^ (f-^0). (3.20) 

For the Wilson loop we get the following result, for both symmetric and 
antisymmetric representations: 


Wt 


N 

00 

1 Q {VX-ir'jRM-l 

k \ 

1 MR 



T k 


IT 


(3.21) 


where Wi is the expectation of the Wilson loop in the fundamental represen¬ 
tation that at strong coupling behaves as [5]: 


Wi = 



^ (\/A-7r)iJM-2 


(3.22) 


This is indeed the correct behavior that we expect at very small k « N. 
The Wilson loop then picks the leading contribution from the term with the 
largest number of traces, due to the large-iV factorization. The result above 
follows from 

tr7^^t/=^ (tr[/)"= +0(tr"-i). (3.23) 

In the high-density regime, symmetric and antisymmetric representations 
behave very differently. For fermions, f grows indehnitely with q and asymp¬ 
totes to f 2 ; 2TT^I‘^\/\n.q at large q. The function h~{f) becomes negative and 
its absolute value grows as the cube of the argument: 


h (f) 2 ; - 


f3 

127r^ 


(f ^ 00 ), 


(3.24) 


which can be seen to match (3.15) in the overlap of the region between the 
two regimes. 
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In the bosonic case, q cannot be bigger than one. When q approaches 
one, the integral in (3.18) converges to a hnite value, indicating that the 
solution for q exists only for sufficiently small f. This corresponds to the 
Bose-Einstein condensation in the analog statistical system. We discuss how 
to compute the symmetric-representation Wilson loop for arbitrary / in the 
next section. 


4 Symmetric representations 

4.1 Bose-Einstein condensation 


The Bose-Einstein condensation for symmetric-representation Wilson loops 
has been discussed for the Gaussian model [23] of AA = 4 SYM |2Sl 1^ - 
as well as in the unitary matrix models [23] • To continue the Wilson loop 
expectation values past the transition point, two prescriptions have been 
used: (i) analytic continuation in / [231129] and (ii) observation that for 
sufficiently large k the rank-/c symmetric representations become equivalent 
to the fc-wrapped fundamental representation [ 2 U] : 

(iiTi+(ti (for sufficiently large A:). (4.1) 


The answer for multiply-wrapped Wilson loops is known exactly at any N 
and k, in the Gaussian model [23]. The scaling limit N ^ oo, k ^ oo can be 
obtained as an approximation to this exact answer [15], and it agrees with 
the symmetric-representation Wilson loop analytically continued far past the 
Bose-Einstein condensation point. The equivalence of the two approaches, 
however, has never been actually proven, and since our model is not Gaussian, 
we prefer to derive all the necessary results from scratch. 

Technically, the condensation happens because the hypergeometric func¬ 


tion in (3.18) has a logarithmic branch point at g = 1: 


1 + K 


2 iq(l, 1 + k; 2 + k; ±q) = - ln(l - g) + 7 + -0(1 + k) + 0(1 - g). (4.2) 


However, this entails no divergences in f, because the logarithmic term inte¬ 


grates to zero in (3.18) and f approaches a hnite limiting value 

00 ^ 2 n-h 


fc = f 1^=1 = TT' 


'E 

n=l 


n 2 e 


-2n 


n 


|2 


= 11 . 


(4.3) 
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c 


V 


Figure 6: The contour deformation in (2.11). 


If f exceeds fc, solution to the saddle-point equation ceases to exist. 

This happens because the chemical potential of a non-interacting Bose 
gas cannot be positive, which translates into u being bigger than /i in onr 
case. The edge behavior of the density of states, p{x) - y/p - x, guarantees 


the convergence of the integral in (2.16) at z/ = /i. Following the analogy with 
statistical mechanics, we expect that the excess density f-fc will condense in 
the gronnd state, identified in onr case with the largest eigenvalne max {aj} = 
a. An important difference to the textbook treatment is the randomness of 
the energy levels which themselves are integration variables with the measure 
defined by the eigenvalne integral ( |2.2[ ). 

To take into account the condensate, we need to deform the contour of 
integration in (2.11) as shown in fig. picking the pole at z/ = a: 




. kLa 


n 


, L{aj—a) 


+ 


(4.4) 


The remaining contonr integral is exponentially small compared to the pole 
contribntion, and can be omitted at large N. The average over the bnlk 
eigenvalues can then be replaced by convolution with the saddle-point den¬ 
sity, while the integral over the largest eigenvalue should be kept as it is: 




CXD /X 

J da P(a) exp kLa-N J dxp{x)\n(^l - 


(4.5) 
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Here -P(a) is the probability to find the largest eigenvalue at point a outside 
the eigenvalue interval: 


F(a) = e 


a^+NLTo 


-2^1 -loop ( ® ® j )) 


(4.6) 


where 


^0 = 


2L/i2 


j dxp{x) In Si_ioop(a -x). 


(4.7) 


The constant Fo is the free energy of the eigenvalue sitting right at the 
edge of the bulk distribution. The probability to find an eigenvalue at a is 
determined by the energy cost of taking an eigenvalue from the edge and 
moving it to a, and it is given by the difference in free energies. 

The probability to pull out one eigenvalue outside the bulk distribution 
is exponentially small, but in the integral (4.5) the negative exponent is 
counteracted by the positive exponent in the thermodynamic weight. The 
largest contribution to the integral hence comes from the saddle point: 


\nW,^ = NLJ^ (a)+ 0{N^). 


This is similar to (2.14), but now the free energy is given by 


F(a) = Fq + /a + — dxp{x) In ^ —_ ^2 


-2^i-ioop(® 2T 

r 


][ _ g L{x-a) 


(4.8) 


(4.9) 


The position of the saddle point is determined by the following equation: 


f = dxp{x) 


S(a - x) - 


, L(a-x) _ 


(4.10) 


with the kernel: 


S{x) E - + - K{x + M) + - K(x - M) - K(x), 
^2 2 


(4.11) 


where we used the explicit form on the one-loop measure (2.4), and defined 

H'{x) 


K{x) = 


H{x) ■ 


(4.12) 
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The first two terms in (4.10) were absent in (2.16) and can be interpreted as 


the density of the condensate. 

Another way to calculate the Wilson loop for f > fc is to continue the 
free energy analytically from f < f,,. Indeed, the logarithmic singularity of 


(4.2) cancels upon integration in (3.18) and f(q'), as a function of q, has 


no singularities aX q = 1. We prove in appendix that the result of the 
analytic continuation gives the same results as the calculation above under 
fairly general conditions, applicable in particular to the model under consid¬ 
eration. The Bose-Einstein condensation consequently does not lead to any 
singularities in the expectation value of the Wilson loop. 

In the universal regime L(a-/i) » 1, the thermal contribution to the free 


energy (the last term in (4.10)) becomes exponentially suppressed and can 
be dropped: 

1 r 2L 

T{a) ^ dxp{x)\nZi_ioop(a-x)-—a^ (4.13) 

A* 

4L 2 r 

f ^ (4-14) 


In this approximation, the /c-symmetric and /c-wrapped Wilson loops have 
the same expectation value: 




(t 


j-QkL^o) 


(4.15) 


because the right-hand side is calculated by equation (4.5) without the second 
term in the exponent 


This result contradicts intuition based on the 
large-factorization, as the Wilson loop now picks the largest contribution 
from the term with the smallest number of traces. 


4.2 Bulk regime 

When a - fi ^ fi, the density is the Wigner distribution and the kernel S{x) 
can be approximated by M‘^B?lx |lj. The answer is then obtained from the 
Wilson loop in the Gaussian model [151123] by simple rescaling: 

NL'^M'^ / _ \ 

^— (k'/I + + arcsinh/tj , (4.16) 
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where 


(4.17) 


K = 


Xf nVXf 


iLfi 2ML 

It is interesting to notice that the consistency of this regime requires the rank 
of representation to scale linearly with the size of the contour. There should 
be a simple explanation to this fact in the dual supergravity picture. 


4.3 Endpoint regime 

We now consider the case when the saddle point a approaches the endpoint 
of the eigenvalue distribution. As before, we introduce the dimensionless 
variables 


V = 


a - /i 

IF 


u = 


yU - X 

IF 


(4.18) 


The density near the endpoint has the form (3.5). It cannot be just substi¬ 


tuted into (4.14) because the integral would then diverge. The strategy is to 


separate the bulk and endpoint contributions to / as follow: 

f = foc+ fep, 


(4.19) 


where foo is computed with the Wigner density and with S{x) approximated 
by M^R^jx: 


/oo X ^ 


M‘^L 

FF 


r , / N 1 4L -T 4LM _ ^ 

/ dxp^(x) - = —\/a^-lR^ - 3 -\F. (4.20) 

J a-x X 0rA4 


Notice that the rank of representation scales linearly with the length of the 
contour as in the bulk regime, but its coupling dependence is different: / ~ 
y-3/4 instead of We thus introduce the rescaled variable 


^ ALM 


such that 


/oo = Vv. 

For the genuine endpoint contribution we hnd: 


(4.21) 

(4.22) 


r / M2i?2 \ 

/ dx[p^{x) -p(x)S'(a - x).). (4.23) 

J \ a-x } 
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In spite of proximity to the endpoint, the argument of the kernel is still very 
big: a-x = M(u+v) » 1, which justihes the use of an approximate expression 

S{a-x) ^ MR^S{u + v) (4.24) 


with 


S{w) = ^ (ta + 1) ln(t(; + 1)^ + ^ ~ 1) - 1)^ - tc Inta^. 


(4.25) 


Writing the endpoint density (2.21), (3.5) as 


p{u) = 


I2M 


TT V h 


F(u), 


(4.26) 


and introducing the function 

oo 

X(v)-j 

0 

we hnd that 


111 \ u + v 


fep = Xiv). 

The function X{v) can be computed by noticing that 

1 1 2 
+ 


tc + 1 w -1 w 


(4.27) 


(4.28) 


(4.29) 


and differentiating (4.27) twice, after which elementary integrations yield: 

6(v - 1) 


X"{v) = 


2Vn 4 


1 _3 

+ - n 2 . 


We thus get: 


X(v) = - ^6'(n-l)(v-l)2 - Vv, 


(4.30) 


(4.31) 


where the integration constants are set to zero in order to get the correct 
behavior at v cx). Interestingly, the function X(n) has a branch-point sin¬ 
gularity outside the eigenvalue support, again due to the resonance appearing 
in the kernel S{x). Combining this result with (4.22), we hnd: 




(4.32) 
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An alternative derivation based on the Wiener-Hopf method [5] is given in 
the appendix [Bj 

Recalling that the equation for / was obtained by differentiating the free 
energy with respect to a, we can get the free energy by integration: 


-^-h/ 


ALM^ 

8LM2 

- 


f u y 

fv-J dw f(w) 




(4.33) 


As before, we would like to express the free energy in terms of /. We need 
then to obtain the inverse function v{f). For n < 1 or equivalently / < 2/3, 
the answer is very simple: 


V = 





ov^r A4 



(4.34) 


At the critical point 



the free energy experiences a phase transition. It is easy to show that T is 
continuous together with its two hrst derivatives, while the third derivative 
diverges at the critical point: 


T = regular + const 



(4.35) 


The phase transition is consequently of the third order. 


4.4 Non-universal regime 

The general scaling dependence of the variable / and of the free energy on 


MR and A have already been discussed in sec. |3]3 

/ - 

R-fif - M{MRy^'^ (4.36) 

Explicit integral representation for the free energy can be obtained by plug¬ 
ging the scaling form of the density (2.22) into (4.9), (4.10). The Fourier 
representation of the the resulting integrals is given in appendix 


23 











5 Conclusions 


The leading-order strong-coupling solution of the SYM* matrix model is 
given by the same Wigner distribution as in the AA = 4 SYM, up to rescaling 
of the coupling A Wilson loop expectation values in the bulk 

regime, as a consequence, are obtained from the Gaussian results by simple 
rescaling. This is quite surprising, as the dual geometry of SYM* is rather 
different from AdS^xS^ away from the boundary. It would be very interesting 
to verify these basic matrix-model predictions by hnding classical D-branes 
solutions in the Pilch-Warner geometry. 

Even more interesting phenomena occur for parametrically small repre¬ 
sentations. The leading order then is just 

(5.1) 

the result expected form the large-Y factorization. It is the next, \/X sup¬ 
pressed term that displays a non-trivial behavior, with phase transitions in 
the representation parameter / = k/N, appropriately rescaled. In the dual 
picture, / maps to the density of electro-magnetic flux on the D-brane world- 
volume. Perhaps the non-analytic behavior of the Wilson loops can be in¬ 
terpreted as phase transitions in the effective held theory on the D-brane 
worldvolume with the electric or magnetic hux playing the role of tunable 
parameter. These phase transitions and quantum phase transitions in the 
matrix model at hnite but large A have the same origin. 

We have mainly concentrated on the inhnite-coupling regime of the the¬ 
ory. The discontinuities in the higher-rank Wilson loops should persist at 
hnite coupling as soon as there are resonance peaks in the densitjj^ which 
is the case for A > A^^^ 35.4 [9]. It would be interesting to completely map 

the phase diagram of Wilson loop expectation values in the A - / plane. 
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A Analytic continuation 

Gonsider an eigenvalue model defined by the partition function 

Z= (A.l) 

i<j 

The eigenvalue density for this model satisfies the saddle-point equation 

^ dy R{x - y)p{y) = V'{x) (A.2) 

with 

(A.3) 

We assume that 

2 

R{x) = - + regular, (A.4) 

X 

such that the kernel in the saddle-point equation is of the Hilbert type, and 
the density has the usual square-root asymptotics near the end-points. It 
will be important for the subsequent derivation that the residue of the kernel 
is exactly two. If it is different from two, the analytic continuation will not 
describe Bose-Einstein condensation correctly. 

The fc-symmetric Wilson loop in this model is determined by the free 
energjj^ 

fj- 

= dxp{x)\n{l - , (A.5) 

The chemical potential u is determined by minimization of the free energy: 

(A.6) 

®Here we set L to one to simplify the notations. The dependence on L can be easily 
recovered by dimensional analysis. 
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which implicitly determines ^ as a function of /. These formulas are literally 
valid for sufficiently small / < fc, such that (A.6) has a real solution. Our 
goal is to analytically continue the free energy past the critical point. 

We first consider the analytic structure of as a function of u. As such, 
it has a cut from -p to /i with the discontinuity 


Disc.F'(x) = 2i7ip(x). (A.7) 

Consequently, near v = p the function has a double expansion that 

combines analytic and square-root terms: 


= / - /c + /i(z^ - /i) + • ■ • + [fi'o + • (A.8) 


The free energy is obtained by integration: 


T{u) = c+{f - p) +-fi{v-pf + ... 




2 2 

- P) + - pf + ■ 


(A.9) 


The constant go must be positive, due to positivity of the eigenvalue density. 
The solution to (A.6) thus only exists for / < fc. 

When fc - f is small, (A.6) can be solved iteratively with the help of the 
expansion (A.8). Substituting the solution back to we get to the first 


approximation: 


nf) = c- 


Uc-ff 


(A.IO) 


While Tiy) has a square-root singularity, the function T{f) is analytic and 
can be continued past fc, which is just an inflection point of the free energy. 
The analytic continuation in / is equivalent to continuing in u through 
the square-root branch cut that changes the sign of fc - f to negative. We 
denote the result of such analytic continuation by 

The function Tiv) has the same continuous part as but the opposite 
discontinuity: 


ConiJ^{x)=!F{x), Disc .F(a:) = -Disc .F(x), (A.11) 

and is completely characterized by these two conditions. 
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We can construct by introducing an auxiliary function, the gener¬ 

alized resolvent: 

W{u) = V{u)-J dyp(y)\YiZ{u-y)-J^o, (A.12) 

-M 

where To is just a constant, to be determined later. The saddle-point equa¬ 
tion (A.2) is equivalent to the two conditions: 

Cont iy'(a:) = 0, Disc ly'(x) = 4i7rp(x), (A.13) 


where the second equation is a consequence of (A.4). This implies, in view 
of (A.7), that 


ContlT(x) = 0, DisclT(x) = 2Disc J’(x). 


(A. 14) 


Strictly speaking, the hrst equation in (A.13) only implies that the continu¬ 
ous part of W{x) is a constant independent of x, but this constant can be 
adjusted to zero by setting 


= r(fi) - J dyp(y) \nZ[y-y) 


(A. 15) 


-/4 


in (A.12). 


It follows from the equations above that 

T{iy)=T{iy)-W{iy). 


(A.16) 


Indeed, this function satishes the right discontinuity condition (A.11). Ex¬ 
plicitly, 

Z(i/ - x) 


T^p) = To +fv-V^p) + J dx p{x) 


In ■ 


(A. 17) 


The function W{v) can be interpreted as the energy cost of pulling one 
eigenvalue out of the bulk of the distribution and placing it at point v. Com¬ 
paring (A.17) with (4.9) we see that the analytic continuation of the free 
energy gives the same result as Bose-Einstein condensation on the largest 
eigenvalue. The Bose-Einstein condensation thus does not lead to any ther¬ 
modynamic singularities in the expectation values of the fc-symmetric Wilson 
loops. 
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B Wiener-Hopf method 

In this appendix we derive the general formulas for the symmetric Wilson loop 
(4.9), (4.10) from the Wiener-Hopf method developed to solve the saddle- 
point equations for the matrix model in jB], assuming that a is very close to 
/i. This covers both endpoint regime and the non-universal regime. As in |5] 
we use the units in which i? = 1 in this appendix, and introduce the endpoint 
variables ^ = fi- x and rj = - a < 0. 

Parametrizing the endpoint distribution as: 


23/2 

P(^) = 


vrp,' 


we can identify three contributions to / in (4.10): 


/ = /oo + /ep + /th, 


(B.l) 


(B.2) 


where 


„ 87r r dx , . 1 + i 

J oo — “T n / Poo\Xj 

A J TT a- 

is the bulk contribution. 


1 + M2 _ 2(M2 + 1) 

X 


TTjJt,^ 




3 

7r/r2 


23/2 


is the end-point contribution and 


/th - 


23/2 fF(Od^ 




I 


3/2 J — 1 


(B.3) 

(B.4) 

(B,6) 


is the thermodynamic term. Here we concentrate on the endpoint part, 
dehned in terms of the function 


oo 

j d(, 




-F(OSU-v) 


(B.6) 


Our main observation is that X(r]) is the same function that was intro¬ 
duced in [B] to solve the integral equation of the matrix model, where it plays 
the role of the remainder function in the Wiener-Hopf method. This function 
















is different from zero only for rj <0 (its Fourier image in analytic in the lower 
half-plane). In the Fourier space, X{u) can be read off from the results in 
[5] (we keep the original notation): 


X(a.) = - 






(B.7) 


E 


(- 1 )’ 


2tt 

/ in<l>{M) 

e M+i 1 


n=l 


nn\ 


/ (M-i)n t 
M+i ) 


/ (M+i) 


e M-i F 


/ {M+i)n\ 
M-i ) 


h + t5)r(-ife) (-‘^ + IS)r(jfc)‘ 


where 


4>{M) = 2Marctan(M) - log (M^ + l) . 


(B.8) 


Let us discuss the analytical structure of this function. The sum has 
two branches of simple poles in the upper half complex plane of u, which 
comes from the G-(cj) dehned in |5]. The prefactor, which is essentially 
1/G_(a;), contains double poles, from F(^) . There is also a square root 
cut in the positive imaginary axis. Therefore, inverse Fourier transforming 
back to the coordinate space for this exact expression is very difficult. As we 
are interested to probe the phase structures at the decompactihcation limit, 
we will study this function at this regime, where the treatment simplihes 
much. 


B.l Universal regime 

To study Tj = 0{M), we deal with u = At the decompactihcation 

limit M -+■ oo, the sum part of (B.7) reduces to (equation (4.6) in [S]): 




sum = 


M^uj 2Msin(^) 
and a naive large M approximation simplihes the prefactor to: 

(-l)h47i-3/2](^3g-5*Ma; (^) 


prefactor = - 
Hence, the full expression is: 

V(u)) = 




uj\/uj - le 


(-l)l/4^3/2JVj-2 / ^ e 

2u}\/oj - ie y 




m 


Muj 

2 


(B.9) 


(B.IO) 


(B.ll) 


29 























We see that X{oj) does not have poles anymore. This is because the simple 
poles got canceled by the zeros of the prefactor. The double poles are not 
relevant at this regime, as T is approximated by its expansion at zero. 
Only the cut remains in the upper half complex plane. 

Straightforward contour integration gives us the expression in the coor¬ 
dinate space, which is: 


X(r]) = [2(-M - - r]) + + 2r])6{-r])^ (B.12) 

3 


Upon rescaling, this is the same as eq. (4.31) in the main text. 


B.2 Non-universal regime 


The solution (B.12) for X{r]) obtained in the regime 0{M) actually approx¬ 


imates well this regime too. 

In Fourier space, following the simplihcations explained in the appendix 
C of [S], X{u) reduces to: 




zu 


2uVuj - ie \ dvr^M Wtt 




(B,13) 


Due to its complicated analytical structure, especially the double poles and 
the logarithmic branch cuts in the exponential, we abstain from computing 
its exact expression in coordinate space. Instead, we study its asymptotic 
behaviors, for a; » 1 and uj « 1. For the former limit, the hrst term of the 


sum in (B.13) dominates, which is the same as the hrst term in (B.ll) from 


the 0{M) regime. For ca « 1, the leading contribution comes from the second 


term of the sum in (B.13), and we expect matching with the regime 0(M), 
which it does. Here are the asymptotic results in the coordinate space: 


\ri\ « 1 : W(u) = -M'^nrj 

|t;|»l; X{v) = 


(B.14) 

(B.15) 


Numerically, we also see that 
large M. 


(B.12) approximates this regime very well at 
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